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of  breaking  to  depth  of  breaking  was  calculated  and  found  to  be 
0.86  for  a composite  of  the  waves  measured;  the  mean  ratio  values 
for  collapsing,  plunging  and  spilling  were  0.84,  0.87  and  0.90,  re- 
spectively. Probability  distributions  were  plotted  for  wave  heights 
and  horizontal  velocities  and  qualitatively  compared  with  Rayleigh 
distributions;  the  wave  distributions  fit  well  but  velocity  did  not. 
The  spectra  of  wave  profile  and  horizontal  velocity  were  calculated 
and  indicated  a narrow  banded  data  set.  The  coherence  values  between 
horizontal  water  particle  velocities  and  wave  profile  were  generally 
high,  indicating  that  the  horizontal  particle  velocities  measured 
were  highly  wave-induced.  The  phase  shift  at  the  peak  energy  fre- 
quency was  about  40°  which  suggested  the  presence  of  reflected  waves. 
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Figure  1.  Figure  of  Variables  for  a Wave  at  the  Breaker 
Position  (Galvin,  1968) . 
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I . INTRODUCTION 


The  surf  zone  is  an  area  bounded  on  the  seaward  side  by  the  point 
where  waves  first  begin  to  break,  and  on  the  landward  side  by  the  point 
of  maximum  run-up  on  the  beach  slope.  There  are  several  classifications 
of  breaking  waves  depending  on  the  wave  steepness,  water  depth,  beach 
slope  and  wavelength.  The  most  generally  accepted  classification  of 
breaking  waves  was  proposed  by  Wiegel  (1954) , and  modified  by  Galvin 
(1968) . In  the  Wiegel  classification,  waves  are  designated  as  spilling, 
plunging  and  surging.  Galvin  (1968)  found  ample  evidence  to  support  a 
fourth  class  of  waves  te:  ed  "collapsing"  (the  type  that  was  most  often 

observed  during  the  period  that  data  were  being  gathered  for  this 
thesis) . 

Spilling  waves  usually  occur  over  gently  sloping  bottoms  and  appear 
once  a wave  crest  becomes  sufficiently  unstable  to  curl  over  slightly 
at  the  peak,  with  the  foamy  portion  "pouring"  down  the  front  slope  of 
the  advancing  wave  form.  Plunging  waves  occur  most  often  on  steeper 
sloping  bottoms.  This  class  of  breaker  is  very  asymmetrical  in  form 
with  the  top  falling  all  the  way  down  to  the  base  of  the  wave  forming 
a large  sheet  splash  and  subsequent  borelike  wave  front.  The  collapsing 
breaker  generally  occurs  on  more  steeply  sloping  beaches  than  spilling 
breakers.  The  collapsing  type  is  very  similar  to  the  plunging  except 
that  it  does  not  break  from  the  highest  point  of  the  wave  but  rather  the 
breaking  portion  becomes  separated  from  the  face  of  the  advancing  form 
at  a position  somewhat  below  the  peak  then  falls  over  in  a fashion 
similar  to  the  plunging  type.  Another  notable  difference  between  the 
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plunging  and  the  collapsing  types  is  that  collapsing  breakers  occur 
at  the  toe  of  the  beach  even  in  the  case  of  very  large  waves,  such  as 
those  observed  at  Carmel  River  Beach. 

The  last  classification  is  the  surging  breaker  which  is  commonly 
observed  on  very  steeply  sloping  beaches.  The  surging  breaker  builds 
to  a peak  but  instead  of  breaking  in  the  usual  way  it  seems  to  subside 
under  its  own  weight  with  the  forward  divergent  momentum  forcing  the 
water  up  the  beach  slope  in  a quiet,  sudden  "surge";  hence  the  name. 

Many  investigators  have  made  attempts  to  define,  derive  or  empiri- 
cally determine  various  breaking  criteria.  Among  those  who  have  pre- 
sented theoretical  developments  are  Boussinesq  (1871) , Lord  Rayleigh 
(1876),  McCowan  (1900),  Rackham  (1952),  Davies  (1952),  Yamada  (1957), 
Lenau  (1966),  Laitone  (1963),  Chappelear  (1959)  and  Dean  (1968).  Of 
the  several  breaking  criteria  proposed  there  are  three  main  categories; 
kinematic,  dynamic  and  geometrical.  A brief  description  of  the  various 
criteria  is  given  in  Table  I. 

There  have  not  been  many  precise  measurements  made  of  breaking 
waves  and  even  fewer  have  been  actually  made  in  the  field.  Most 
measurements  have  been  made  in  laboratory  wave  tanks. 

For  this  thesis  breaking  wave  profiles  and  the  particle  velocities 
beneath  them  were  measured  in  the  field  on  a very  steeply  sloping 
(0.11)  beach.  The  waves  varied  in  height  from  approximately  0.22  to 
0.92  meters. 

The  objective  was  to  measure  probability  distributions  of  and 
n^,  for  breaking  waves  in  the  surf  zone  and  the  water  particle  velocities 
beneath  the  waves.  Measurements  were  made  on  a steep  beach  to  ascertain 
the  shape  dependence  on  breaking  wave  parameters.  The  ultimate  objective 
of  this  and  related  studies  is  to  devise  a scheme  for  predicting  the 
wave  height  and  water  depth  at  breaking. 
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TABLE  I.  THEORETICAL  BREAKING  CRITERIA 


KINEMATIC 


DYNAMIC 


GEOMETRICAL 


V ^ Cfa,  i.e.,  the  limiting  value  for  the  water 
particle  velocity  at  the  wave  crest  is  greater 
than  the  wave  phase  speed.  If  the  particle 
velocity  exceeds  this  value,  it  would  separate 
itself  from  the  waveform. 


1)  Maximum  vertical  particle  acceleration  can- 
not exceed  the  gravitational  acceleration  (p 
constant,  surface  pressure  constant) . The  water 
particle  leaves  the  surface  vertically  for 
larger  accelerations. 

2)  Limiting  pressure  gradient  is  zero. 


The  surface  slope  of  the  wave  steepens  as  the 
wave  shoals.  The  limiting  value  for  this  slope 
is  infinite  which  is  a vertical  face.  Beyond 
this,  the  water  particles  fall  forward  ahead 
of  the  wave. 
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II.  MEASUREMENTS 


A.  EXPERIMENT  SITE 

During  the  period  of  3 May  to  12  May  1976  an  experiment  was  conducted 
at  Carmel  River  Beach  in  the  Monterey,  California  area  (see  Figure  2). 
Four  days  were  spent  putting  the  equipment  in  place.  Because  of  large 
wave  conditions  the  task  was  both  arduous  and  time  consuming.  Removing 
the  equipment  was  much  easier  due  to  the  fact  that  a hurricane  off  Baja 
California  produced  a large  swell  and  associated  breakers,  one  of  which 
carried  away  most  of  the  equipment,  thus  terminating  the  experiment. 

The  Carmel  River  Beach  site  was  chosen  for  several  reasons.  First, 
it  is  located  in  an  embaynent  which  causes  the  waves  entering  the  bay  to 
be  directionally  filtered  and  refracted  by  the  geometry  of  the  bay  so 
that  narrow-banded  waves  impinge  almost  perpendicular  to  the  beach. 
Therefore,  in  one  respect,  the  description  of  the  waves  is  greatly  sim- 
plified. 

The  second  reason  for  choosing  the  Carmel  River  Beach  site  was  to 
study  waves  breaking  on  a steep  beach  slope.  The  resulting  wave  system 
was  not  at  all  simple  because  of  the  presence  of  strongly  reflected 
waves.  The  steep  slope  produced  breakers  primarily  of  the  collapsing 
type  with  lesser  numbers  of  spilling,  plunging  and  surging  breakers. 

Figure  3 shows  the  beach  profile  and  instrument  location  at  Carmel 
River  Beach  on  10  May  1976.  Table  II  gives  the  characteristics  of  the 
beach. 
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Distance  in  Meters 

Figure  3 . Beach  Profile  and  Instrument  Location  at  Carmel  River  Beach 

10  May  1976 


TABLE  II . BEACH  AND  WAVE  CHARACTERISTICS  ON  10  MAY  1976 


Location: 


Carmel  River  Beach 


Beach  Slope: 

Sand  Composition: 

Sand  Median  Diameter: 


Breaker  Type: 
Wave  Period: 


Wave  Height: 


9:1 


Quartz  feldspar 

0.6  mm  (after  Bascom,  1946) 


Collapsing  and  Plunging 
10.56  seconds 


Up  to  3.5  meters 
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B.  INSTRUMENTS 


The  wave  gauges  were  of  the  capacitance  type  constructed  of  3/8  inch 
outside  diameter  stainless  steel  rod  covered  with  1/16  inch  wall  thick- 
ness polypropylene  tubing  that  fit  tightly  over  the  rod;  these  wave 
gauges  proved  to  be  sufficiently  sturdy  to  withstand  the  severe  force  of 
the  surf  zone.  The  only  real  problem  that  occurred  with  the  gauges  was 
one  of  water  proofing  the  electrical  contacts.  The  gauges  operate  on 
the  principle  that  a change  in  capacitance  is  noted  as  a change  in  plate 
dimension  occurs.  In  these  gauges  the  3/8  inch  steel  rod  and  sea  water 
act  as  the  two  plates  and  the  polypropylene  cover  acts  as  the  dielectric. 
As  the  water  level  rises  around  the  gauge,  as  in  a passing  wave,  there 
is  a change  in  the  capacitance  which  is  sensed  by  a self-contained 
battery  powered  transistorized  circuit  that  produces  a linear  change  in 
the  output  voltage.  The  circuit  was  designed  by  McGoldrick  (1969)  and 
housed  in  a watertight  container  mounted  on  the  support  tower. 

A static  calibration  of  the  wave  gauge  was  accomplished  in  the  labora- 
tory by  lowering  the  gauge  into  a deep  tank.  Using  the  battery  powered 
circuit  and  short  connecting  leads  to  reduce  induced  capacitance,  an 
accuracy  of  + 0.005  meters  was  established.  A plot  of  voltage  output  as 
a function  of  depth  of  emersion  was  plotted  and  the  slope  of  the  curve 
was  used  to  determine  the  conversion  constant  from  volts  to  meters.  The 
calibration  plot  is  shown  in  Appendix  A. 

The  flow  meter  used  was  an  EPCO  Model  6130  Electromagnetic  Ocean  Cur- 
rent Meter  which  works  on  Faraday's  principle  of  electromagnetic  induc- 
tion. The  meter  measures  water  velocity  in  two  orthogonal  directions. 

The  flow  meter  was  dynamically  calibrated  by  oscillating  it  in  a water 
tank.  This  procedure  was  used  because  the  characteristics  were  expected 
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to  be  different  under  rapidl(y  changing  flow  conditions  than  those  pro- 
duced by  steady  flow  (see  Appendix  B) . The  problems  encountered  with 
the  flow  meters  were  caused  by  the  extremely  harsh  conditions  to  which 
they  were  subjected  such  as  abrasion  by  waterborne  sand,  gravel  and 
various  other  forms  of  debris,  including  the  experimenters  themselves 
when  knocked  about  by  the  waves. 

Two  Interstate  Electronics  Corporation  Model  SDP  201  pressure  trans- 
ducers were  placed  in  5.5  meters  of  water  at  a distance  of  47.5  meters 
directly  seaward  of  the  outermost  tower.  The  pressure  sensor  is  a 
small,  unbonded  strain  gauge  bridge.  Direct  current  excitation  for  the 
bridge  is  supplied  by  a voltage  regulator  located  in  the  transducer 
housing.  The  sea  pressure  is  coupled  by  a neoprene  diaphragm  to  a sili- 
cone fluid  filling  the  interior.  One  part  of  the  transducer  is  exposed 
directly  to  the  interior  fluid  and  the  other  is  connected  to  a chamber 
which  is  connected  to  the  interior  tluid  by  a length  of  capillary  tubing. 
This  arrangement  acts  as  a hydraulic  filter  developing  a reference 
pressure  which  is  an  average  value  of  the  external  sea  pres. sure.  By  the 
action  of  this  filter  the  transducer  senses  rapid  pressure  fluctuations 
only  and  slow  changes  such  as  tides  are  lost  through  the  hydraulic  filter. 
The  pressure  transducers  were  put  in  place  by  divers,  and  the  signals 
were  cabled  to  the  beach. 

The  instruments,  with  the  exception  of  the  pressure  transducers, 
were  attached  to  a 3.6  meter  high  tower  constructed  of  steel  pipe  with 
an  oustide  diameter  of  6.3  cm.  A 1.0  meter  diameter  base  plate  was  placed 
about  0.6  meters  from  the  bottom  of  the  tower.  This  configuration  allowed 
the  towers  to  be  sunk  approximately  0.65  meters  into  the  sand  and  provide 
firm  support.  The  towers  were  also  supported  by  steel  guy  wires  fastened 
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to  blade  anchors  driven  into  the  sand.  At  low  tide  the  first  tower  was 
placed  as  far  as  possible  seaward.  The  seaward  extent  was  limited  only 
by  the  ability  of  the  experimenters  to  wade  around  in  1 meter  of  water 
and  sledge  an  anchor  into  the  bottom  while  waves  of  1 to  1.5  meters 
broke  over  their  heads.  However,  after  considerable  effort  the  seaward 
tower  was  placed.  The  second  tower  was  located  shoreward  of  the  first 
in  a line  perpendicular  to  it. 

The  capacitance  wave  gauge,  along  with  a stiffening  rod  which  also 
acted  as  the  ground  lead,  was  mounted  on  the  outer  tower.  A specially 
designed  bracket  was  used  to  hold  the  wave  gauge  20  cm  away  from  the  tower 
and  thus  avoid  any  possible  interference  between  the  tower  and  the  gauge. 
The  flow  meter  was  mounted  directly  below  the  wave  gauge  and  close  to 
the  bottom. 

The  instruments  were  tied  to  the  shore  facilities  by  coaxial  cables 
run  over  the  tops  of  the  towers  to  a support  tower  on  the  beach  berm. 
Located  at  the  support  tower  was  a push  button  switch  that  actuated  a 
one  volt  signal  generator.  The  switch  was  installed  so  that  an  operator 
using  a prearranged  code  could  "mark"  each  wave  that  broke  directly  on 
the  wave  gauge  according  to  visually  estimated  breaker  type  (spilling, 
plunging,  collapsing  or  surging) . The  mark  in  the  form  of  an  electrical 
impulse  was  recorded  and  used  to  categorize  breaker  type.  All  signals 
were  recorded  using  a Vidar  Corporation  32-channel  digital  data  acquisi- 
tion system.  An  8-channel  Brush  strip  chart  recorder  was  used  for 
continuous  monitoring  of  the  instruments. 
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III.  ANALYSIS 


The  data  collected  at  Carmel  River  Beach  on  10  May  1976  were 
selected  for  analysis.  Table  II  contains  the  conditions  prevailing 
during  the  experiment.  The  data  include  a record  of  instantaneous 
water  level,  horizontal  water  particle  velocities  and  identified  and 
classified  waves  that  broke  directly  on  the  wave  gauge.  The  digitized 
wave  records  were  converted  to  analog  records  and  plotted  using  an 
IBM  360  computer  in  conjunction  with  a CALCOMP  model  765  plotter. 

A mean  water  level  was  computed  using  a running  mean  technique.  Seg- 
ments of  the  total  record,  each  384  seconds  long,  were  plotted.  The 
reason  for  this  particular  record  length  was  to  put  a maximum  number 
of  waves  on  a minimum  number  of  plots  and  still  preserve  visual  resolu- 
tion. Horizontal  water  particle  velocities  underneath  the  waves  were 
plotted  in  a similar  fashion. 

The  plots  were  manually  analyzed  in  the  following  way.  First  the 
marked  breaking  waves  were  measured  for  height  at  breaking.  Next  the 
value  of  depth  at  breaking  was  determined  and  referenced  to  mean  water 
level,  which  Galvin  (1968)  proposed  as  more  reasonable  than  being 
referenced  to  still  water  level.  The  ratio  of  the  breaker  height  to 
the  depth  at  breaking  was  computed  as  was  the  ratio  of  horizontal 
water  particle  velocity  to  calculated  wave  speed  at  breaking.  Prob- 
ability distributions  were  calculated  for  each  of  the  parameters  and 
compared  to  a Rayleigh  distribution.  The  variances,  standard  devia- 
tions and  average  values  were  calculated  for  each  of  the  parameters. 
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Spectra,  coherence  and  phase  between  waves  and  horizontal  veloci- 
ties were  also  calculated.  The  data  were  sampled  every  0.00375  seconds 
resulting  in  a Nyquist  frequency  of  2.0833  Hz.  Three  30-minute  records 
were  analyzed.  Using  a maximum  lag  time  of  5%  of  the  record,  the  band- 
width resolution  is  0.00728  Hz.  Computing  the  spectra  in  this  manner 
results  in  40  degrees  of  freedom  for  each  spectral  estimate.  The  80% 
confidence  limits  for  40  degrees  of  freedom  using  a chi-square  distri- 
bution are  between  0.73  and  1.30  of  the  measured  power  spectral 


estimates . 


IV.  RESULTS 


A.  SPECTRA 

Analysis  of  the  field  data  (Appendix  C)  collected  at  Carmel  River 
Beach  on  10  May  1976  showed  that  of  the  398  waves  recorded,  59  broke  on 
the  gauges  and  these  waves  ranged  in  size  from  0.22  meters  to  0.92 
meters.  The  point  of  breakingjOr  depth  at  which  a wave  breaks, is  func- 
tionally related  to  the  wave  height.  Since  wave  height  can  be  considered 
as  a random  variable  it  follows  that  the  point  of  breaking  will  tend  to 
be  random  as  well.  Because  the  point  of  breaking  is  random, only  about 
15%  of  the  waves  on  10  May  1976  broke  directly  on  the  measurement  towers. 
The  larger  waves  that  were  present  were  breaking  "outside"  or  seaward  of 
the  gauge  and  those  which  were  of  smaller  size  were  breaking  "inside” 
or  shoreward  of  the  gauge. 

Figure  4 is  a plot  of  the  log  spectral  density  of  the  waves  and  the 
horizontal  water  particle  velocity.  The  plot  represents  the  data  from 
one  30  minute  segment  of  a 90  minute  data  set  recorded  at  Carmel  River 
Beach  on  10  May  1976.  The  spectrum  is  narrow  banded  with  the  peak  energy 
at  0.0947  Hz,  corresponding  to  a 10.56  second  period.  There  is  a second 
peak  in  the  spectrum  approximately  half  the  primary  frequency  and  repre- 
sents an  approximate  sub-harmonic  of  0.0947  Hz.  Guzza  and  Davis  (1974) 
state  that  in  a theoretical  analysis  of  edge  waves  excited  by  incoming 
waves  the  prominent  edge  wave  mode  is  the  first  sub-harmonic  of  the  pri- 
mary frequency.  Finally,  the  spectral  value  at  about  0.012  Hz  probably 
represents  surf  beat  which  generally  falls  in  the  range  from  0.01333  Hz 
to  0.0083  Hz  or  periods  of  75  to  120  seconds. 
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PHASE  COHERENCE  LOG  SPECTRAL.  DENS  I TV  N**Z-SE'C 


Figure  4.  Log  Spectral  Density  for  Waves  and  Horizontal 
Velocities,  Phase  and  Coherence. 
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The  coherence  values  between  waves  and  horizontal  velocity  were 
moderately  high,  greater  than  0.70  indicating  the  velocities  are  pri- 
marily wave  induced  over  a frequency  range  of  0.0510  to  0.2622  Hz. 

But,  for  the  rest  of  the  frequency  range  the  values  are  quite  small. 

The  crests  of  breaking  waves  tend  to  lead  the  horizontal  particle 
velocities  in  the  body  of  the  wave.  There  was  approximately  a 40  de- 
gree phase  lag  between  the  wave  crest  and  the  horizontal  velocity  at 
the  peak  frequency.  The  phase  difference  between  waves  and  velocities 
would  be  expected  to  be  zero  for  progressive  waves.  The  large  value 
is  presumably  due  to  the  presence  of  reflected  waves. 

B.  BREAKING  WAVES 

The  values  of  the  height  at  breaking  were  measured  directly  from 
the  plots  of  the  wave  record.  Wave  record  plots  produced  by  the 
CALCOMP  675  were  first  matched  with  plots  of  the  wave  marker  (Figure  5) 
and  each  wave  of  interest,  that  is  one  which  broke  directly  on  the  tower, 
was  marked.  The  vertical  distance  between  the  maximum  and  minimum  water 
surface  elevations  about  a zero  up-crossing  at.  the  breaker  point  was 
measured  and  this  value  used  as  the  breaker  height.  This  definition  of 
breaker  height  was  in  most  instances  adequate;  however,  in  some  cases 
a problem  arose  when  trying  to  analyze  the  record  using  this  definition. 
The  problem  manifested  itself  whenever  there  were  waves  in  the  record 
that  were  obviously  composed  of  harmonics  or  secondary  waves,  as  in  the 
case  of  wave  33  in  Figure  5.  Values  of  the  parameters  for  the  total 
record  are  presented  in  Appendix  C. 

It  was  found  that  the  mean  value  of  the  amplitude  of  the  waves  at 
breaking  was  a little  more  than  two- thirds  of  the  wave  height  at  breaking, 
or  nK  = 0.68  H.  . If  one  considers  these  values  relative  to  the  mean 
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water  level  it  is  seen  that  the  height  of  the  wave  crest  above  mean 
water  level  is,  in  the  mean,  1.5  times  the  depression  of  the  trough 
below  mean  water  level.  The  mean  value  of  1.5  for  the  crest  to  trough 
ratio  results  from  averaging  over  the  whole  record;  some  values  were 
much  less  but  in  some  extreme  cases  the  ratio  was  as  great  as  5.5. 

There  was  not  enough  data  to  conclusively  state  that  the  ratios  were 
much  greater  for  one  particular  wave  type  than  another  although  visual 
observations  suggested  that  the  crest  to  trough  ratio  was  the  greatest 
for  collapsing  breakers. 

The  value  of  the  depth  of  breaking  (h^)  was  taken  to  be  the  depth  of 

the  water  from  the  mean  water  level  to  the  bottom.  This  value  was 

measured  at  the  point  of  zero  up-crossing  following  the  trough  of  a 

wave  immediately  preceding  one  of  the  marked  waves.  The  mean  water  level 

was  chosen  as  the  reference  level  in  determining  the  value  of  h,  because 

b 

it  gives  a better  representation  of  the  actual  water  level  at  any  given 
instant  than  still  water  level.  The  average  water  level  was  calculated 
using  a running  mean  techique  averaged  over  800  data  points,  or  24 
seconds . 

Even  the  mean  water  level  is  not  always  a particularly  good  reference 
for  the  value  of  the  height  of  breaking.  It  was  observed  that  some  of 
the  waves  actually  broke  at  a point  where  the  water  depth  in  the  backwash 
was  almost  zero. 

Waves  shoal  rapidly  on  a steeply  sloping  beach  face.  The  first  and 

most  obvious  occurrence  is  the  steepening  of  the  wave  from  an  almost 

« 

indistinguishable  crest  line  in  deeper  water  to  a very  definite  crest  and 
trough  system  closer  to  shore.  The  wave  celerity  slows  and  the  leading 
edge  at  the  base  of  the  trough  may  even  slow  to  a stop  relative  to  the 


28 


backwash  running  down  the  beach  slope.  This  supports  the  idea  that  the 
swash  can  be  a dominant  factor  on  steep  beaches.  At  the  stationary 
"foot"  of  the  wave,  drawdown  and  gravity  are  working  together  to  supply 
water  to  the  growing  "wall"  or  waveface  (see  Figure  6) ; but  this  is 
not  happening  fast  enough  and  therefore  a water  deficit  develops.  Tur- 
bulence and  the  resulting  trapped  air  is  carried  to  the  base  of  the 
building  wave  and  up  the  face  creating  voids  below  the  wave  crest.  Ver- 
tically moving  water  particles  also  gain  forward  momentum  from  the 
advancing  crest,  causing  them  to  separate  from  the  wave  face  before 
reaching  the  crest.  Thus  the  wave,  lacking  support  and  having  forward 
motion  in  the  upper  regions,  gives  the  appearance  of  a plunging  breaker 
that  has  had  its  support  removed.  As  Galvin  (1968)  pointed  out  these 
breaking  waves  occur  at  the  toe  of  the  beach,  or  the  lowest  reach  of  the 
backwash.  The  toe  of  the  beach  is  the  only  place  that  this  can  occur 
since  a water  deficit  is  necessary  to  have  a collapsing  breaker. 

This  phenomenon  does  not  show  up  on  wave  records  because  a mean 
reference  or  still  reference  is  used,  which  by  its  very  nature,  masks 
short  term  extremes.  Since  the  collapsing  mechanism  occurs  in  the  trough 
region  immediately  preceding  a breaking  wave,  perhaps  the  instantaneous 
water  level  below  the  trough  should  be  thought  of  as  the  depth  of  breaking. 


C.  DEEP  WATER  WAVES 

A mean  value  of  the  initial  deep  water  wave  height  (FT)  was  calculated 
from  the  spectrum  of  the  offshore  pressure  record.  The  pressure  record 
was  converted  to  a surface  elevation  using  linear  wave  theory  for  which 
pressure  at  depth,  AP,  due  to  the  waves  is  given  by 
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The  wave  profile  is 


ri  = a cos  (kx  - at) 


(2) 


so  that 


AP  - (pg  cosh  ^(^z)) 
cosh  kh 


(3) 


The  term  in  the  brackets  expresses  the  spectral  transfer  function  H^(f) 
relating  the  surface  profile  spectrum  to  the  pressure  spectrum 

S (f)  = |h  (f)  |2  S (f)  (4) 

1 1 

The  deep  water  wave  amplitude,  aQ,  can  be  related  to  the  amplitude 
at  the  measurement  depth,  a^.  Assuming  conservation  of  energy  and 
applying  linear  wave  theory 
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where  kq  = initial  deep  water  shoaling  coefficient 

= shoaling  coefficient  at  point  of  breaking. 

The  term  in  the  brackets  is  the  spectral  transfer  function, 
relating  the  deep  water  spectrum  to  the  surface  profile  spectrum  at 
any  depth. 


(7) 
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but  the  mean  height  at  breaking  from  the  data  was  found  to  be  0.42 
meters,  which  was  less  at  breaking  than  in  deep  water.  From  theory  this 
should  not  be  the  case.  As  a wave  shoals,  thq  wavelength  shortens  and 
the  height  increases  while  the  period  remains  the  same.  It  was  first 
postulated  that  reformed  waves  were  being  measured  resulting  in  the  low 
breaking  wave  heights  observed,  but  this  idea  was  rejected.  It  was 
decided  that  reflected  waves  very  likely  were  the  reason  for  the  apparent 
contradiction . 

Reflected  waves  would  be  expected  to  be  important  since  the  slope 
at  Carmel  River  Beach  at  the  time  of  observation  was  a very  steep  0.11. 
The  data  further  indicated  that  the  horizontal  shoreward  velocities  are 
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about  the  same  or  somewhat  less  than  the  seaward  swash  velocities.  The 


difference  in  these  velocities  could  be  an  effective  mechanism  in  the 
generation  of  reflected  waves.  For  a reflected  wave  off  an  idealized 
plane  boundary  the  nodal  point  is  a distance  of  L/4  from  the  reflecting 
surface  (see  Figure  7) . The  approximate  distance  from  where  the  mean 

I 

water  level  intersects  the  beach  to  the  measurement  tower,  x,  can  be 
determined  from  the  beach  geometry,  such  that 

tan  6 = - (11) 

x 

where  tan  B = beach  slope  and  h = mean  water  depth  = 0.49. 

Solving,  x = 4.5  meters. 

The  wavelength  can  be  approximated  by 


L = CT  ~ /gn  T . (12) 

The  average  period  was  10.6  seconds  so  that 

L/4  5.8  meters  . 

The  calculations  show  that  the  measurement  point  was  near  the  node  of  the 
reflected  wave  system  which  would  result  in  relatively  small  wave 
heights  at  the  measurement  tower. 

D.  BREAKING  WAVE  PARAMETERS 

Waves  tend  to  arrive  at  the  beach  in  groups  which  surfers  call  "sets". 
The  groups  of  relatively  large  waves  periodically  supply  more  water,  or 
momentum  flux,  to  the  beach  and  cause  a super-elevation  or  increased 
set-up  inside  the  surf  zone  while  the  group  persists.  As  the  last  waves 
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of  the  group  break,  the  water  level  begins  to  fall  and  remains  character- 
istically low  until  the  next  "set"  arrives.  This  periodic  rise  and  fall 
of  the  surface  elevation  is  referred  to  as  "surf  beat". 

A rise  in  the  mean  water  level  was  observed  to  occur  shortly  after 
the  arrival  of  a "set"  or  group  of  large  waves.  The  amount  of  deviation 
of  the  depth  at  breaking  calculated  about  the  mean  value  of  the  depth  at 
breaking,  ranged  from  minus  10.25%  to  plus  12.28%.  This  variation  in 
water  level  is  one  explanation  why  waves  of  such  varied  sizes  could  break 
at  essentially  the  same  point  in  space. 

The  ratio  of  breaker  height  to  depth  was  calculated.  This  ratio  has 
long  been  used  as  a criterion  for  wave  breaking.  In  other  studies  the 
ratio  of  H^/h^  has  been  used  over  a wide  range  of  bottom  and  slope  con- 
ditions in  both  laboratory  and  field  measurements.  The  theoretical  values 
range  from  about  0.73  (Boussinesq,  1871)  to  1.03  (Packham,  1952).  It  is 
generally  found  in  field  work  that  these  are  reasonable  values. 

In  1944  field  data  were  gathered  on  breaking  waves  at  the  Scripps 
Institution  of  Oceanography  pier  (Munk,  1947) . This  data  constitutes 
a major  portion  of  the  relatively  few  records  of  breaking  waves  actually 
measured  in  the  field.  The  values  of  H^/h^  ranged  from  0.67  to  0.82. 

In  this  thesis  the  mean  value  of  the  ratio  H^/h^  was  found  to  be  0.86 
which  is  within  the  range  of  theoretical  values.  However,  this  value 
is  larger  than  field  measurements  made  by  Scripps  and  smaller  than  the 
0.96  reported  by  Woods  Hole  (Gaughan,  et  al.,  1973).  The  SIO  data  along 
with  data  taken  on  10  May  1976  at  Carmel  River  Beach  is  given  in  Figure  8. 
Woods  Hole  data  are  not  included  because  no  slope  information  was  avail- 
able. Figure  8 is  a plot  of  H^  versus 

the  SIO  report  and  from  this  experiment.  The  values  of  wave  height  in 


h^  and  contains  data  from  both 
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the  SIO  data  were  much  larger  than  those  in  this  report.  For  reference 
purposes  the  theoretical  linear  slope  of  0.78  from  solitary  wave  theory 
is  plotted  in  the  figure. 

A comparison  between  the  field  data  obtained  by  SIO  and  this  experi- 
ment reveals  that  many  individual  values  of  H, /h  for  the  Carmel  River 

b b 

Beach  were  considerably  larger  in  magnitude.  The  beach  slope  during  the 

SIO  experiments  was  0.04  and  less  while  the  slope  at  Carmel  River  Beach 

was  0.11.  The  difference  in  beach  slope  presumably  accounts  for  the 

larger  height  to  depth  ratios  observed. 

The  ratio  of  H./L,  was  also  calculated.  The  wavelength  at  breakinq, 
b b 

L^,  was  determined  from  the  expression  relating  wavelength  to  celerity 
and  period. 


'b  b 


(13) 


where  wave  speed  is  approximated  as, 


C 


b 


(14) 


This  is  a rough  approximation  to  the  value  of  L^,  but  from  previous  work 
done  by  Galvin  (1975) (where  he  used  two  wave  towers  to  measure  earlier 
at  this  beach)  it  was  found  that  the  above  formula  slightly  overestimated 
the  celerity,  again  presumably  due  to  a reflected  wave  system. 

The  classification  of  breaker  types  according  to  Battjes  is  shown  by 
the  curves  in  Figure  9.  Van  Dorn's  (1977)  data  shown  in  the  figure  is 
according  to  Battjes  (1974)  reinterpretation  of  Galvin's  (1968)  observa- 
tions. The  value  of  H./L,  found  in  the  present  study  is  also  plotted  in 
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Figure  9,  and  indicates  that  the  field  evaluation  of  the  breaker  type 
agrees  with  the  Battjes  classification. 

E.  PROBABILITY  DISTRIBUTIONS 

Probability  distributions  were  plotted  for  Hfa  and  p as  shown  in 
Figures  10  and  11,  respectively.  In  each  case  these  distributions  at 
least  qualitatively  resemble  a Rayleigh  distribution.  Barber  (1130) 
was  the  first  to  relate  a Rayleigh  distribution  to  the  heights  of  sea 
waves.  Longuet-Higgins  (1955)  showed  theoretically  that  the  distribu- 
tion of  wave  height,  H,  for  the  case  of  a narrow  banded  spectra  is 
given  by  the  Rayleigh  distribution, 


P (H)  = -JL-  e-(H/2)2/202 
(20)  2 


(15) 


Here  we  apply  the  Rayleigh  distribution  to  breaking  waves  by  letting  H 

2 

be  the  wave  height  at  breaking;  o is  the  sample  variance  of  the  surface 
profile . 

Based  on  data  available  at  the  time  (mostly  from  offshore  pressure 
sensor  records)  Longuet-Higgins  (1952)  found  a difference  between 
measured  and  theoretical  values  for  a Rayleigh  distribution  of  only  8% 
or  less.  The  strict  assumptions  made  in  deriving  the  theoretical  prob- 
ability distribution  makes  this  agreement  surprising  and  it  may  indicate 
that  the  narrowness  of  the  spectrum  is  not  as  important  as  assumed. 

Since  the  work  by  Longuet-Higgins,  Koexe  and  JeBruyn  (1964)  and 
Goodknight  and  Russel  (1963)  have  f ,u*_  for  deep  water  waves  the 

goodness  of  fit  value  of  the  Rayleigh  i. s tr ibution  is  very  high.  Battjes 
(1974)  also  points  out  that  the  old  dictum  proposed  by  Rayleigh  himself. 
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Figure  10.  Probability  Distribution  of  Wave  Heights  and 

Corresponding  Theoretical  Rayleigh  Distribution. 
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Field  Data 
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Figure  11.  Probability  Distnbut.on  of  Wave  Amplitudes  and  Correspond- 
ing Theoretical  Rayleigh  Distribution. 


"The  basic  law  of  the  seaway  is  the  apparent  lack  of  any  law,"  seems 
somewhat  in  doubt  based  on  these  findings. 

The  following  explanation  is  given  by  Battjes  (1974)  as  the  reason 
why  wave  heights  are  to  a good  approximation  Rayleigh  distributed  even 
for  waves  that  are  non-linear  such  as  at  breaking  and  do  not  possess  a 
narrow  spectrum.  Even  though  breaking  waves  have  higher  crests  and 
shallower  troughs,  the  increase  in  crest  compensates  for  the  shallowing 
of  the  trough  and  H values  remain  essentially  unchanged.  Taking  the 
definition  for  crest  elevation  as  being  the  greatest  elevation  between 
successive  zero-crossings,  it  was  assumed  that  the  crest  elevation  cor- 
relates strongly  with  the  immediately  following  largest  trough  depth 
between  successive  zero-crossings.  Battjes  observed  that  Jahns  and 
Wheller  (1972)  reject  this  assumption.  They  point  out  that  for  a wide 
spectrum  a relatively  high  crest  is  on  the  average  followed  by  a trough 
that  is  not  relatively  deep.  They  argue  that  application  of  a Rayleigh 
distribution  would  lead  to  an  overestimation  of  the  probability  of  ex- 
ceeding some  arbitrary  value  for  the  relatively  high  waves  if  the 
spectrum  were  not  narrow. 

It  was  found  in  the  Carmel  River  Beach  data  that  the  Rayleigh  distri- 
bution generally  tends  to  underestimate  the  probability  of  exceeding  some 
arbitrary  value  of  nb  and  overestimates  the  probability  of  exceedence 
for  values  of  H^.  This  result  tends  to  agree  with  Jahns  and  Wheller's 
argument. 

The  variance  used  in  calculating  the  Rayleigh  distribution  was  de- 
rived from  the  total  record  of  instantaneous  water  level  whereas  the 
breaking  wave  height  data  only  includes  those  waves  that  broke  on  the 
tower.  Hence,  the  calculated  variance  leads  to  an  cverestimation  of  the 
probability  of  the  wave  height  exceeding  some  arbitrary  value. 
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Probability  distributions  were  also  calculated  for  horizontal  water 
particle  velocities  measured  under  the  breaking  waves,  and  are  shown  in 
Figure  12.  Distributions  were  calculated  for  motion  in  both  directions, 
i.e.,  seaward  and  shoreward.  These  were  compared  to  the  Rayleigh  dis- 
tribution, but  it  was  found  that  the  probability  of  exceeding  some 
arbitrary  value  is  considerably  underestimated. 

The  distributions  of  the  velocities  in  the  seaward  direction  are 
strongly  negatively  skewed  showing  more  large  maximum  velocities  than 
small  values.  The  shoreward  velocities  are  negatively  skewed  but  not 
to  the  same  degree  as  the  seaward.  The  larger  seaward  velocities  would 
be  expected  on  a steep  beach  because  the  seaward  velocity  is  the  sum 
of  backwash  due  to  the  force  of  gravity  plus  the  particle  velocity  in 
the  trough  of  the  wave.  The  Rayleigh  distribution  is  positively  skewed 
and  therefore  does  not  represent  the  distribution  of  the  observed 
velocities. 

( + ) 

The  probability  distributions  for  the  ratio  of  V /C,  for  both 

max  b 

collapsing  and  plunging  breakers,  shown  in  Figure  13,  indicate  that  waves 

break  over  the  range  of  values  from  0.5  to  1.5.  The  kinematic  criteria 

states  that  breaking  will  occur  when  V > C,  (Table  I).  The  breakers 

<+>  ~ b 

that  occurred  for  values  of  V /C,  < 1.0  indicate  that  other  processes 

max  b 

may  control  breaking,  such  as  the  water  deficit  mentioned  earlier. 

Figure  14  is  a composite  of  all  the  Vmax/Cb  ratios  computed  on  10  May 
1976. 
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Rayleigh 

(-) 

V 


Figure  12.  Probability  Distribution  of  Shoreward  Horizontal  Velocity 

under  a Breaking  Wave,  Seaward  Swash  Velocity  and  Correspond- 
ing Theoretical  Rayleigh  Distribution. 


V.  CONCLUSIONS 


The  following  conclusions  are  derived  from  the  field  measurements 
made  on  Carmel  River  Beach  on  10  May  1976.  The  predominate  breaker 
types  were  collapsing  and  plunging.  The  beach  slope  was  0.11  indicating 
a steep  slope. 

1.  Very  little  difference  between  the  ratio  of  H^/h^  for  the 
various  breaker  types  was  measured. 

2.  The  mean  ratio  of  H^/h^  calculated  from  the  data  was  0.86  and 
is  close  to  generally  accepted  theoretical  values. 

3.  Breaking  wave  distributions  of  and  are  qualitatively  de- 
scribed by  the  Rayleigh  distribution.  The  probability  of 

exceedance  is  overestimated  in  the  case  of  H,  and  underesti- 

b 

mated  in  the  case  of  n.  • 

b 

4.  Distribution  of  the  onshore  and  offshore  water  particle  veloci- 
ties in  breaking  waves  is  negatively  skewed  and  therefore  not 
described  by  the  Rayleigh  distribution 

5.  The  mean  value  of  the  wave  height  at  breaking  may,  in  the  case 

of  reflected  waves,  be  less  than  the  initial  deep  water  wave 

height.  This  will  only  occur  close  to  a nodal  point. 

(+) 

6.  The  value  of  the  ratio  of  V /C,  was  computed  for  all  the  waves 

max  b e 

that  broke  on  the  tower  and  in  many  cases  the  ratio  was  less  than 

1.0.  This  was  most  prominent  for  the  collapsing  breakers.  Since 

the  kinematic  breaking  criterion  states  that  V > C,  at  break- 

max  — b 

ing,  backwash  and  wave  reflection  must  be  important  in  describing 
breaking  criterion  on  steep  beaches. 
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7.  The  breaking  mechanism  in  collapsing  breakers  appears  to  be 
the  turbulence  present  in  the  base  of  the  wave  caused  by  a 
water  deficit  in  the  swash  zone. 


APPENDIX  A 


The  calibration  of  the  wave  gauge  along  with  its  associated  lineariz- 
ing circuit  was  carried  out  in  the  NPS  physics  laboratory  deep  acoustics 
tank.  Figure  15  shows  a calibration  plot  for  wave  gauges. 

CALM  is  the  slope  of  the  curve  in  Figure  15  in  meters/volt 

CALA  = height  of  gauge  above  bottom  related  to  the  intercept  of  the 
curve. 

The  Engineering  Physics  Company  water  current  meter  was  calibrated 
by  oscillating  it  in  a water  tank,  as  described  by  Steer  (1972) . Steer 
obtained  a calibration  factor  of  0.610  meters/sec/volt  for  the  probe  used 
in  this  experiment. 
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APPENDIX  B 


This  appendix  contains 
(+)  (-) 

for  Hb/hb,  nb,  V , V , 


a collection  of  probability  distributions 
(-) 

and  V /C  . 
max  b 
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Collapsing  Breakers 
Plunging  Breakers 
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Figure  16.  Probability  Distribution  of  H^/h  for  Collapsing  and 
Plunging  Breakers. 
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Figure  19.  Probability  Distribution  of  for  Collapsing  and 
Plunging  Breakers. 
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Figure  20.  Probability  Distribution  of  V for  Collapsing  and  Plunging  Breakers 
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Figure  22. 
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APPENDIX  C 


This  appendix  consists  of  a table  containing  the  data  collected 

during  the  Camel  River  Beach  experiment  on  10  May  1976.  The  table 

contains  all  the  data  taken  directly  from  the  plots  as  well  as  some 

ratios  computed  from  these  data.  For  these  data,  m = 0.11  and 

H ' =0.49  meters, 
o 
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During  Carmel  River  Beach  (Thornton-Hultstrand) 
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